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We theoretically study the efltects of spin-orbit coupling on spin exchange in a low-density 
Wigner crystal. In addition to the familiar antiferromagnetic Heisenberg exchange, we find gen- 
eral anisotropic interactions in spin space if the exchange paths allowed by the crystal structure 
form loops in real space. In particular, it is shown that the two-electron exchange interaction can 
acquire ferromagnetic character. 
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To a first approximation, the electrons in a Wigner 
crystal are localized in space through their mutual 
Coulomb repulsion. The crystal electrons are thus dis- 
tinguishable through their location on the crystal lattice. 
For any finite interaction strength, however, the electrons 
are able to tunnel through the localizing Coulomb bar- 
rier, loosing this distinguishability. For spinless particles, 
such tunneling is largely inconsequential, since the elec- 
tron configurations before and after tunneling are equiv- 
alent. Because they carry spin, however, a process in 
which two crystal electrons exchange their positions does 
have an effect: it exchanges the spins on the respective 
lattice sites. While the charge configuration of the crystal 
is still nearly static, its spins acquire some dynamics. 

In the absence of spin-orbit interaction (SOI), the cor- 
responding two-electron exchange interaction takes the 
Heisenberg form, JSi • S2, by spin-rotational symme- 
try. Because typically the lowest-energy orbital wave 
functions are symmetric under particle interchange, the 
spin exchange is generically antiferromagnetic (J > 0), 
as a consequence of the Pauli principle. This state- 
ment has been established rigorously for one-dimensional 
(ID) many-electron systems with velocity- and spin- 
independent interactions [1]. The SOI, however, has 
long been known to modify the Heisenberg form of the 
exchange Hamiltonian by, e.g., effectively canting the 
participating spins. The resulting Dzyaloshinsky-Moriya 
(DM) interaction [2] of the form D • (Si x S2), with some 
structure-dependent vector D, was initially studied in the 
context of weak ferromagnetism in materials that other- 
wise were expected to be antiferromagnetic, but has since 
been appreciated in many other contexts. 

This paper continues the recent discussion of the SOI 
in mesoscopic and Wigner-crystal exchange processes 
[31 m |5j |6]. For definiteness, we focus our attention on 
clean strongly interacting quasi-lD systems [Tj. Two- 
electron as well as ring exchange interactions in single- 
and double-row quasi-lD systems without SOI were ex- 
tensively discussed in Refs. [Elin]- We show that the SOI 
qualitatively enriches the electronic behavior in ID by 
causing anisotropics in the exchange Hamiltonian. Our 
calculation is performed using the path-integral instan- 




FIG. 1: (Color online) Electrons confined to move along the 
X axis undergo a geometric spin transformation as they ex- 
change their positions. The shaded area marks the transverse 
confinement of width w. As the electrons tunnel through the 
mutual Coulomb repulsion barrier, they deviate laterally by 
the distance r ~ w{'w / asY^^ , in the center-of-mass refer- 
ence frame shown here. Sketched are the classical trajectories 
following the rim of the inverted potential. If the left elec- 
tron chooses to follow the upper path, the right electron is 
forced to simultaneously go along the lower path and vice 
versa. U± is the SOI SU(2) transformation Q starting at the 
origin of the xy coordinate system, moving to the left, then 
(counter)clockwise along the upper (lower) rim and back. The 
total area enclosed by the two trajectories is A^ r"^ . 

ton picture of particle exchange [HI [10]. In the low-density 
limit, as <C h, where h = is the average inter-electron 
separation and as = eh? /me^ is the Bohr radius, the 
dominant electron paths follow classical trajectories in 
the inverted potential. The SOI is then naturally cap- 
tured by purely geometric SU(2) spin transformations 
along the electronic exchange paths [11]. See Fig. [l] for a 
schematic. We thus explore the non-Abelian Aharonov- 
Casher [12] variant of physics whose Aharonov-Bohm 
counterpart was studied extensively in Refs. [T5] . 

To illustrate our approach and the key findings, we 
consider two spin-1/2 electrons confined to move in the 
xy plane, subject to the two-dimensional (2D) effective- 
mass Hamiltonian (denoting by hats the spin structure): 

H=^ [{-ihVj + Aj)V2m + V,] + C/12 . (1) 
i=i,2 

Here, Vj = {dx^.dy.) and Aj = A(rj) = {A^^Ay) is 
the two-component 2x2 vector potential, which can 
be expanded in the basis of the three Pauli matrices 
{ax: o'j/i o-^) = &. The SU(2) vector potential A describes 
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the SOI. The Pauli equation for electrons in vacuum, for 
example, has A cx Vl^ x o", while in 2D electron systems 
with asymmetric confining potential along the z axis, the 
SOI is usually dominated by the so-called Rashba interac- 
tion of the form: A = az x o-. Further, U\2 = U(y:\ — "c-i) 
is the effective two-particle repulsion and Vj = Vivj) is 
the 2D confinement potential along the x axis. The two- 
electron exchange then becomes a quasi- ID scattering 
problem: In the center-of-mass reference frame, electrons 
are incident from opposite directions and are scattering 
near the origin, Ar = ri — r2 ~ 0. The electron transmis- 
sion through the repulsion potential [/(Ar) in the low- 
density limit is well approximated by WKB tunneling. 

We start by assuming a homogeneous Rashba SOI. In 
the strictly ID limit, when the lateral motion along the 
y direction can be completely neglected, it is possible to 
gauge out the SOI by the unitary transformation 

= 6-*^^="/'*. (2) 

When applied to each electron separately, this gauge 
transformation removes the vector potential A from the 
Hamiltonian. It is, however, not possible to generalize 
this to 2D, since and Ay in general do not commute. 
We, nevertheless, choose to transform our 2D Hamilto- 
nian ([ij as if ^ W^HW. When the electrons are suf- 
ficiently distant from each other and they stay close to 
y = due to the V{y) confinement, the transformation 
([2]) then does remove the SOI from the problem. Only 
during the short time intervals when the electrons tun- 
nel through each other do they deviate appreciably from 
the X axis and undergo an additional SOI-induced spin 
transformation (see Fig. fll. It is this remnant spin pre- 
cession that the gauge ^ allows us to focus on. To 
simplify our discussion, we assume that the spin preces- 
sion length Iso = Trh/a is much longer than the width 
w — y^h/mu! of the wire (assuming quadratic lateral 
confinement V = mtj^y^/2), while the Bohr radius qb 
is much shorter than w: <C w <C /so- In the opposite 
limit, w <^ ub, electrons are confined very tightly in the 
lateral direction, forcing them to tunnel "head on" (9j and 
thus invalidating the simple instanton calculation of the 
geometric SOI phase. In the case of 2D and 3D Wigner 
crystals, it would be sufficient to require that <C b, 
which guarantees that the electrons avoid head-on tun- 
neling. In the case of a quasi- ID model, it is convenient 
(although not essential) to assume 6 » r ~ w{w/aBy^^, 
which allows one to neglect the tendency of ID crystals 
to form zigzag arrangements 

Before proceeding, it is instructive to recall how one 
constructs the exchange Hamiltonian in the low-density 
limit without SOI using the instanton method. The 
strength of the exchange is parametrized by a positive 
real- valued parameter J, which is determined by the Eu- 
clidian action Se along the minimal-action path exchang- 
ing two electrons [TD]: J = f3nuJo^/st^j2^e-'^''/'', 



where ujq is the characteristic attempt frequency (cor- 
responding to the effective electron confinement along 
the wire) and /3 is a numerical prefactor of or- 
der unity. The classical minimal action path corre- 
sponds approximately to the rim of the inverted po- 
tential C/ij(Ar) = -[V{Ay)/2 + U{Ar)] (again assum- 
ing quadratic confinement) in the center-of-mass frame: 

~ j'^i, dly/2m*\UEil)\- J is thus essentially the 
WKB tunneling amplitude for a particle with the re- 
duced mass TO* = to/2 through the potential barrier Ue 
along the classical trajectory parametrized by I. The 
authors of Ref. [9 performed a detailed study of the Eu- 
clidian action in quasi-lD wires, obtaining in the regime 
of as <C w, r <C b: Se/It- ~ TjyJVfoB — k{w [asY^^ , where 
71 and K are numerical prefactors of order unity. The first 
term on the right-hand side is the action for the head-on 
tunneling of electrons and the second term is its reduc- 
tion due to lateral excursions that allow electrons to avoid 
each other, as sketched in Fig. [l] 

In the absence of the SOI, J parametrizes the usual 
antiferromagnetic coupling between the two spins [5]: 

iJAF = J(i «) i + o- (g) 0-) , (3) 

making the convention that the left (right) operator in a 
tensor product L®R acts on the left (right) spin. The op- 
erator X = (l®l + a-®&)/2 simply exchanges the spins: 
-'^Isi, S2) = |s2, Si), where |si, S2) is the spinor wave func- 
tion with si (S2) corresponding to the left (right) spin, 
respectively. In the case of multiple exchange trajecto- 
ries, J stands for the sum of all contributions. We now 
include SOI effects. During the position exchange, each 
electron that is moving along path C undergoes the path- 
dependent SU(2) transformation [IT] 

[/ = Tcexp-''^c''"--A/;i^ 

where Tq is the contour-ordering operator and the inte- 
gral runs over the classical exchange paths closed along 
the X axis, as shown in Fig.[l] We choose the point where 
the two electrons pass each other to be at x = 0, see 
Fig. [1] and integrate (counter)clockwise along the upper 
(lower) classical exchange trajectory. The SOI is assumed 
to be weak in comparison to the interparticle repulsion, 
such that the effect of the spin dynamics on the instanton 
orbital trajectories may be neglected. 

The effective spin Hamiltonian is finally obtained as: 

Hs ^ J{ul ®U- + u^® u+){i ®i + a®a)i2. (5) 

The operator ([/I ® U+)X implements the spin transfor- 
mation along the clockwise exchange trajectory, when the 
left electron deviates from y — Q toward positive y, mov- 
ing to the right, while the right electron is correspond- 
ingly pushed toward negative y, moving to the left. The 
other term u]^ ® [/_ accounts for the counterclockwise ex- 
change path. Hg is verified to be Hermitian through the 
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identity X{U+ (g) Ul) = {Ul_ (g) U+)X. Notice that elec- 
trons whose exchange occurs at x = xq instead of a; = 0, 
as in Fig. [T| acquire an additional rotation of their inter- 
action Hamiltonian by W{xq): Hg W''' {xq)HsW{xq), 
with W applied to both spins. 

We now employ Eq. Q to calculate the spin exchange 
coupling in the presence of the SOI, Eq. ([T|). Parametriz- 
ing U± — u± — iu± ■ (T with real- valued scalars u± and 
vectors u±, such that u^ + u^ — 1, the Hamiltonian ([s]) 
acquires the form (omitting a constant offset): 

Hs/J =(?i+M_ — U+ • U_)o-i • (T2 + d • {&i X 0-2) 

+ 2(a • 6-i)(a • ^2) - 2(b • &i){h ■ &2) , (6) 

where d = ii+u_ + u_u+ parametrizes DM, a = (u_|- + 
u_)/2 Ising antiferromagnetic, and b = (u+ — u_)/2 
Ising ferromagnetic interactions. The operators &i and 
&2 SiCt on the left and right spin, respectively. 

Specializing to the Rashba case and expanding the 
transformation matrices in 7 = ^/A{a/h) (where A is 
the total area of the loop formed by the two classical 
trajectories in Fig.[T]), we have [TT| : 

U±^lTi-/^^z + il^^±-o- + 0{-/'^). (7) 

Since the cubic-order in 7 contribution to U± depends on 
the exact shape of the exchange loop, we parametrized 
it by the dimensionless spin-space vectors v± . The same 
is true also of the higher-order terms. Substituting the 
expansion ^ into Eq. ([6| gives: 

Hr^'/J - 6-1 • ^2 - 7'v • (<3-i X &2) ~ 2j^ai,a2. , (8) 

where v = v+ + v_. In addition to the leading anti- 
ferromagnetic coupling, we find a DM interaction [2] at 
order 7'^ and a ferromagnetic Ising coupling along the 
z axis at 0(7"*). In the approximation ([8|, we retained 
only the leading in 7 terms separately for Heisenberg, 
DM, and Ising interactions. If the wire is mirror sym- 
metric with respect to the xz plane, we have ~ u_ 
and u = — A/u+, where M stands for the mirror image, 
so that u+ + u_ = u+ — MU4. (X y. This means, in 
particular, that d oc y in Eq. ([6| and the DM interaction 
Hbm oc y • ((Ti X 0-2) is of the general form allowed by 
the mirror symmetry, to all orders in 7. Furthermore, a 
Rashba system is mirror symmetric with respect to the 
xy plane (combined with flipping a —a). This con- 
strains the DM coupling to be odd and the Heisenberg 
and Ising terms even in a, in accordance with Eq. ([s]). 

A ferromagnetic Ising coupling of the form (TizfT2z is 
also expected as a consequence of correlated orbital quan- 
tum fluctuations, which produce van der Waals-type spin 
interactions via SOI [6 . A form of the exchange similar 
to our Eq. (|8|, consisting of the Heisenberg, Ising, and 
DM pieces has been reported before in Refs. [SlIlllS]- An 
analogous result was also found for the RKKY interac- 
tion mediated by itinerant electrons in the presence of 



the Rashba SOI [S]. In contrast to our Eq. ([6]), however, 
Refs. [3j m [3 [6] predicted a spin exchange of the form: 

^ J&i • <3-2 + Tn • i&i X <3-2) +r(n ■&i)in- 6-2) , (9) 

parametrized by a single vector n. This stemmed from 
a hidden SU(2) symmetry of the exchange Hamiltonian 
in Refs. [31 [3] (when the Hamiltonian can be written as 
a Heisenberg exchange between canted spins) and from 
a spin-rotational symmetry around n in Ref. Such 
symmetries are not assumed in our calculation based on 
Eq. ([5|, allowing for the more general Hamiltonian ([6|. 

To elucidate the anisotropic spin structure of the ex- 
change Hamiltonian ([6]), we distinguish two effects of 
the SOI. First, the SO coupling cants the participat- 
ing spins through a spin rotation along the exchange 
path. We have already gauged out the main part of 
that rotation by means of our transformation W, but 
deviations of the exchange paths from the x axis add 
another piece, contributing to U±. A canting of two 
spins that participate in the usual Heisenberg exchange 
[31 [51 [5] by a rotation angle 9 around the direction n, 
iJ„/J = {Vi&,v}) ■ {V^&2V2), where Vj = e-™"^^ ''/^^ 
produces an exchange Hamiltonian of the form Q , with 
T — — Jsin 29 and F = J(l— cos 9). Hn, however, has the 
same eigenvalues as the isotropic Heisenberg exchange 
Hamiltonian ^ [51 [Ti]: A mere canting of two spins 
should not be considered a real anisotropy, although the 
resulting Hamiltonian has a DM and an Ising pieces. 

Our exchange geometry has a 2D character with two 
different exchange paths corresponding to transforma- 
tions U±. This results in an exchange Hamiltonian Hg 
whose eigenvalues differ from those of H^p, Eq. ([s]). In 
order to quantify this second and more interesting ef- 
fect of the SOI, and to disentangle it from the effects of 
a mere canting of the participating spins, we bring the 
exchange Hamiltonian into a standard form through 
local spin rotations. To this end, we first rewrite Hg, 
Eq. ([6|, as Hs/J — o'ih&2, in terms of a real-valued 
3x3 tensor h . Through a singular- value decomposition 
(SVD), we then bring h into the form h = RDT^ , with 
rotation matrices R,T £ S0(3) and a real-valued diago- 
nal matrix D. The diagonal matrices D are the desired 
standardized representation of exchange Hamiltonians. 
We find for our Hamiltonian ([6]): D^^ = 1, Dyy = 1, 
and Dzz — u+u^ + u+ • The sign of Dzz deter- 

mines whether the exchange interaction in the rotated 
spin coordinates has antiferromagnetic {D^z > 0) or fer- 
romagnetic (Dzz < 0) character: Rotations by tt around 
the z axis can be used to flip the signs of D^x and Dyy, 
such that all entries of D have the same sign as Dzz- 

Note that our exchange Hamiltonian ^ always has a 
pair of degenerate singular values \Djj\, for any 11+ and 
J7_, realizing an XXZ model in rotated spin coordinates. 
This general statement is easiest understood by express- 
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ing the square of the Hamiltonian ([5|, 
4^2/ j2 ^ [ljllJ_ ® lJ_ljl + jyi_U+® U+Ul) + 2 , (10) 

in a basis where the transformations (j\_U- acting on the 
left spin and U- Ul, acting on the right spin are diagonal. 
The only spin operator occurring in the expression for 
is then aiz<^2z- Evidently, is severely constrained: it 
has an Ising form in proper coordinates. This allows us 
to draw conclusions about Hs itself, after writing it in 
the general form Hs/ J = a'ih&2 + c, where we restored 
the constant c that was omitted in Eq. ([6]). In the spin 
basis where h is diagonal also the tens or representing 
Hs is diagonal and, according to Eq. ( 10 1, only one of its 



elements is nonzero. Straightforward algebra then shows 



that this implies D 
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XX 



so that \Dx 



\Dyy\ (up 



to a permutation between three cartesian axes). 

We are now ready to compare the exchange Hamilto- 
nian Eq. obtained in Refs. O HI El [6] to our 
Hs , Eq. (IgI . Representing Hi^^ by a corresponding tensor 
we find that also h^^^ has at least one pair of iden- 
tical singular values. H^^^ thus realizes an XXZ model in 
a rotated spin basis, just as our Hg, Eq. (jo]). Eq. (jo]) has 
been derived in the strongly interacting limit under the 
assumption that only a single pair of (time-reversed) ex- 
change paths contributes. For the case that the inverted 
repulsion potential has more than one saddle point, one 
finds that the exchange Hamiltonian after a SVD will in 
general take the form of an XYZ model. The same holds 
true if the orbital dynamics in the x direction allows par- 
ticle exchanges at a variety of locations xq. The resulting 
exchange Hamiltonian Hg {xo)HsW{xo) averaged 
over Xq may then also realize a Heisenberg exchange with 
two independent anisotropy axes, an XYZ model. 

Let us turn briefly to electronic systems that are ef- 
fectively ID through the confining potential V{y). In 
the absence of SOI, the generic antiferromagnetic ex- 
change Hamiltonian Haf there results in the expected 
Luttinger- liquid behavior at low energies [15j . Devi- 
ations from Haf, however, can have profound conse- 
quences. In Refs. [16], for example, spin-charge sepa- 
ration described by a new universality class has been 
found for a ID Bose gas with ferromagnetic Heisenberg 
exchange. One may thus expect important implications 
of our Eq. ([b]) for strongly interacting quantum wires. It 
was shown in Refs. [9j that the effects of the coupling of 
the exchanging electron pair to surrounding electrons in 
ID may be systematically included in the instanton ap- 
proach, producing only small corrections in typical lim- 
its. In a uniform wire without SOI, the two-electron 
exchange Hamiltonian then carries over to the many- 
electron system, HJ'"° = Hg. With SOI, however, that 
is no longer true, as the Hamiltonian H^"^ acquires a 
position dependence through our gauge transformation 
|2|. Tjwiro = W^HgW. The pairwise exchange interac- 



tion thus becomes z-dependent (unless Hg is invariant 
under W), which couples the orbital motion to the spin 
dynamics. We furthermore remark that even if W has no 
effect on H^^'^'^, one generally cannot bring the exchange 
Hamiltonians into a diagonal form for all pairs of neigh- 
boring spins simultaneously by a proper choice of spin 
bases. The rotations in the above SVD decomposition 
are in general incompatible for consecutive spin pairs. 

In conclusion, we have analyzed spin exchange in effec- 
tively ID interacting electron systems at low density. In a 
two-electron system, the resulting exchange Hamiltonian 
can be brought into the form of an anisotropic Heisen- 
berg model. The model can have both antiferromagnetic 
and ferromagnetic character. We have discussed general 
conditions for the degree of anisotropy of the resulting 
exchange. Both XXZ and XYZ models may be realized. 
Our results may have profound implications for interact- 
ing many-electron systems in one dimension. 
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